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Dynamic properties of molecular motors whose motion is powered by interactions 
with specific lattice bonds are studied theoretically with the help of discrete-state 
stochastic "burnt-bridge" models. Molecular motors are depicted as random walkers 
that can destroy or rebuild periodically distributed weak connections ("bridges") 
when crossing them, with probabilities pi and p2 correspondingly. Dynamic proper- 
ties, such as velocities and dispersions, are obtained in exact and explicit form for 
arbitrary values of parameters pi and p2- For the unbiased random walker, reversible 
burning of the bridges results in a biased directed motion with a dynamic transition 
observed at very small concentrations of bridges. In the case of backward biased 
molecular motor its backward velocity is reduced and a reversal of the direction of 
motion is observed for some range of parameters. It is also found that the dispersion 
demonstrates a complex, non-monotonic behavior with large fluctuations for some 
set of parameters. Complex dynamics of the system is discussed by analyzing the 
behavior of the molecular motors near burned bridges. 



INTRODUCTION 



In recent years an increased attention has been devoted to investigations of molecular 
motors, also known as motor proteins, that are crucial in many cellular processes They 
transform chemical energy into the mechanical motion in non-equilibrium conditions. For 
most of molecular motors their motion along linear molecular tracks is fueled by the hydroly- 
sis of adenosine triphosphate (ATP) or related compounds. It was suggested that a different 
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mechanism is employed to power the motion of a protein collagenase along collagen fibrils 

a a. 

It probably utilizes the collagen proteolysis, cleaving the filament at specific sites. 
As the collagenase molecule is unable to cross the already broken bond, it leads to the bi- 
ased diffusion along the filament. However, full understanding of mechanisms of collagenase 
motion is still not available. 

It was proposed that a good description of the collagenase dynamics could be provided 
by the so-called "burnt-bridge model" (BBM) 3, S 3, B, S E, B Q. In this model, the 
motor protein is depicted as a random walker that translocates along the one-dimensional 
lattice that consists of strong and weak bonds. While the strong bonds remain unaffected if 
crossed by the walker in any direction, the weak ones (termed "bridges") might be broken 
(or "burnt") with a probability < < 1 when crossed in the specific direction, and the 
walker cannot cross the burnt bridges again, unless they are restored, which can occur with 
probability < ^2 ^ 1- In Refs. 6|, l7| an analytical approach was developed which per- 
mitted us to derive the explicit formulas for molecular motor velocity V{c,pi) and diffusion 
constant D{c,pi) for the entire ranges of burning probability < pi < 1 and concentration 
of the bridges < c < 1 which were also confirmed by extensive Monte Carlo computer 
simulations. This theoretical method has been applied to several problems with periodic 

n Q 

bridge distribution. However, the results in |6J, |7| have been obtained only for irreversible 
bridge burning (bridge recovery probability was taken to be P2 = 0), and also for unbiased 
random walker between bridges (equal forward and backward transition rates). In present 
work, we generalize our approach to allow for the possibility of bridge recovery as well as 
unequal hopping rates on the sites between bridges. It is more realistic to consider systems 
with reversible action of motor proteins since they are catalysts that equally accelerate both 
forward and backward biochemical transitions [1]. 



MODEL 



According to our model, we view a motor protein as a random walker moving along 
an infinite one- dimensional lattice with forward and backward transition rates being u and 
w correspondingly, as illustrated in Fig. 1. The lattice spacing size is set to be equal to 
one. The lattice is composed of strong and weak bonds. There is no interaction between 
the random walker and strong bonds, however crossing the bridge in the forward direction 
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Figure 1: A schematic picture of the motion of a molecular motor in the reversible burnt-bridge 
model. Thick solid lines depict strong links, while thin solid lines represent periodically distributed 
weak links (bridges). Dotted lines are for already burnt bridges. 

(from left to right) leads to its burning with the probability pi, while the particle moves 
with the rate u. After the weak link is destroyed, the walker is assumed to be on the right 
side of it. When the particle is trying to cross a broken bond, the bridge can be recovered 
with the probability p2i while the particle moves to the left with rate w. It is assumed that 
initially, at t = 0, all bridges are intact. 

The details of breaking weak bonds in BBM have a strong effect on the dynamic properties 
of motor proteins Q]. There are two different possibilities of bridge burning. In the first 
variant (the so called "forward BBM"), the weak bond is broken when crossed from left to 
right, but the intact bridge is not affected when the particle moves from right to left. Thus 
the bridge recovery may occur if the walker attempts to cross a burnt bridge from right to 
left. In the second variant (named "forward-backward BBM"), the weak link is destroyed if 
crossed in either direction 0, [sl. Both variants are identical for pi = 1, however for pi < 1 
the dynamics is different in two burning scenarios, as was shown in p2 = case JgI], although 
mechanisms are still the same. For reasons of simplicity, below we will only consider forward 
BBM, even though forward-backward BBM can also be solved using the same method. 

There are five parameters that specify the dynamics of molecular motors in BBM: the 
probabilities pi, p2, the concentration of bridges c, as well as transition rates u and w. 
The dynamic properties of the walker are also strongly influenced by the distribution of 



weak bonds js]. Below we will study the case of periodically distributed bridges, when their 
concentration is c = and the weak bonds are located between the lattice sites with the 
coordinates kN — 1 and /cA^with integer k (see Fig. 1). This description is more realistic 
for coUagenases' dynamics [2, 13|. The model below will be studied using continuous time 
analysis as it better describes chemical transitions in motor proteins 
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Figure 2: Reduced kinetic scheme for continuous-time forward BBM with bridge recovery (only 
transition rates not equal to u, w are shown). The origin is the right end of the last burnt bridge. 
Parameters are described in detail in the text. 

DYNAMIC PROPERTIES FOR BBM WITH BRIDGE RECOVERY 



Velocity 



n 



To find the walker's velocity, we generalize the method used in |6i] (for irreversible bridge 
burning, i.e. p2 = 0) to allow for non-zero probability p2- We introduce a probability Rj{t) 
that the random walker is found j sites apart from the last burnt bridge at time t. The 
probabilities Rj{t) arise if the system is viewed in moving coordinate frame with the last 
burnt bridge always at the origin, as illustrated by a reduced chemical kinetic scheme shown 
in Fig. 2. 

The dynamics of the system is determined by a set of master equations: 
dRkN+iit) 



dt 



uRkN+i~l{t) + wRkN+i+i{t) - (m + w)RkN+i{t), 



(1) 



for /c = 0, 1, 2, ■ ■ ■ and i = 1, 2, ■ ■ ■ , - 2; and 

dRkN+N-l(t) 



dt 



uRkN+N-2{t) + Wp2f{k + l)Ro{t) + wR[k+l)N{t) -{U + w)RkN+N-l{t), (2) 



dR^ 



{k+l)N 



it) 



dt 



[1 - Pi)uRkN+N-l{t) + wRi^k+l)N+l{t) - (m + w)R(k+l)N{t), 



(3) 
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with k = 0,1,2, ■■ ■ for both Eqs. ([2]) and ([3]). Also at the origin we have 
dR (t) °° 

= PiuJ^i^kN-iit)] + wRiit) - uRoit) - wp2Ro{t). (4) 

k=l 

In Eq. ([2]) we introduced a function f{k) as a probabihty that next to the last burnt bridge is 

oo 

k periods to the left from the last burnt bridge. It satisfies the condition ^ f{k) = 1, which 

fc=i 

is reflected in Eq. (jl]). The system of equations is to be solved in the stationary-state 

limit (at large times) when dRj{t)/dt = is satisfied, and we denote Rj{t — oo) = Rj in 
what follows. By definition, it can be argued that 

rf,^ PlRkN-1 ,r\ 

^ PlRkN-1 
k=l 

Solving the system can be facilitated by rewriting Eq. ([5]) in a more convenient form. 

To this end, we note that based on the results from Refs. 0, 8| it is reasonable to assume 
that RkN+i is of the form 

RkN+^ = y'WiT), (6) 

where y and W are some functions of pi, p2, u, w and N. Furthermore, y is i and k- 
independent, while W depends on i (but not on k). The advantage of the ansatz ([6]) is that 
it leads to a simpler form of f{k): 

fik) = y'-\l-y), (7) 

as follows from Eq. ([5]). We proceed to solve Eqs. ([I])-(11D with f{k) in Eq. ([2]) given by the 
expression ([7]). 

Introducing a parameter (3 = u/w, it can be easily verified that 

RkN+^ = Cl{k)+C2{k)P' (8) 

with arbitrary Ci{k), C2{k) solves Eq. ([1]). Utilizing Eq. ([8]), Ci{k), C2{k) can be expressed 
in terms of functions RkN and RkN+i (first two points of the period). Then Eq. leads to 

w 

RkN+i = RkN 1 ~ Z^* '^k, (9) 

U — W ^ 

where we defined 

Afc = RkN+l — RkN- (10) 
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In expression ([9]) it was assumed that k = 0, 1, 2, ■ ■ ■ and i = 0,l,---,A^ — 1. Parameters 
RkN and are to be determined from Eqs. ([2]) and ([3]). Substituting Eq. ([9]) in ([2]) and 
(ED yields 



uR 



kN 



UW 



U — W 



[1 - Afc + wp2y''{l - y)Ro + wR^k+m -{u + w)R, 



+ + nj)^^ [1 - A, = 0, (11) 



u — w 



mil 

(1 - pi)uRkN - (1 - Pi) [1 - P""-'] Afc + wR 



u — w 



(k+l)N 



W 



U — W 



[1 - (3] Afc+i - (u + w)R(k+i)N = 0. 



Eq. ( ITTi) also yields 



u — w r-. t/ 



w(l-/3^) 
Substituting ([I3D into Eq. (USD results in 



[i?fc7v - Rik+i)N - P2y (1 - y)Ro]- 



aRkN + bR(^k+i)N + ci?(fc+2)7v + id + ey + fy'^)Roy^ = 



where 



a = (1 



1-/3 



N-l 



1-P 



N 



= —u + (1 — PljU— — w 



c = w 



d = p2{l-pi)it 



1-/3^ 

1-/3 
1-/5^' 

1 - /j^-^ 



1-P 



AT' 



and 



1 - (3^-^ 1 - 13 

e = -P2(l - Pl)u + Wp^Y^: 

\-(3 



Using Eq. ([6]), it follows that 



and Eq. flMj) turns into 



RkN = Roy^, 

{c + f)y^ + ib + e)y + ia + d) = 0. 



(12) 



(13) 



(14) 

(15) 
(16) 
(17) 
(18) 
(19) 

(20) 

(21) 

(22) 
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Therefore 



-{b + e) - v/(& + e)2-4(a + d)(c + /) 



(23) 



2(c+/) 

where we selected the solution of ( !22|) such that < y < 1. With the help of ( !2T1) we obtain 
from Eq. ([I3D 



A, 



u — w 



Roy\l-y){l-p2). 



Thus ([9]) yields 



kN+i 



Roy' 



l-(l-p2)(l-2/) 



(24) 



(25) 



Eq. fl25|) with y given by fl23|) solves the system of Eqs. in the stationary state limit. 

We note that although the Eq. (jlj) was not used to find this solution, it was numerically 
verified that every equation in the system ([I]) -(jlj) is indeed solved by expressions (l25l) and 

m. 



Parameter Rq needed to find the velocity is found from Eq. (!25j) combined with the 

oo N-1 

normalization condition Yl Yl RkN+i = 1, which produces 

k=0 i=0 



i?0 = (l-l/)|iV-(l-P2)(l-l/) 

The mean velocity of the walker is given by 



N 



1 



1-/3^ 1-p 



(26) 



V = ^(wj - Wj)Rj = {u- wp2)Ro + 

j=0 



oo N-l 



^ ^(m - w)RkN+i 



,k=0 i=Q 



(u — w)Ri 



0, 



(27) 



which results in a simple relation, 



V = w{l- P2)Ro + (u- w). 



(28) 



In Eq. ( I28i) . Rq is given by ( l26l) with y from the expression (J23j). 

It can be shown that in the limit of m — > 1, w — * 1, and P2 Eq. fj28l) reproduces 
the result obtained earlier in the Ref. jf| for the BBM with u = w = 1 and p2 = 0. Also, 
Eq. (j28j) simplifies considerably in the limiting case of pi = 1 (deterministic bridge burning) 
when y = 0. In the case of pi = 1 we obtained V{u,w,p2,N) in [oj] using the Derrida's 
method lOj] and our general result given in Eq. ( l28l) agrees with it in the pi ^ 1 limit, as 
was numerically verified. 
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Diffusion Coefficient 

The diffusion coefficient is found by generalizing the method developed in [t] (where we 
found dynamic properties of the random walker in BBM with ^ = ^ = 1,^2 = and 
periodic bridge distribution), allowing for < p2 < 1 and u ^ w. We define PkN+i,mit) 
as the probability that at time t the random walker is located at point x = kN + i {i = 
0, 1, ■ ■ ■ , — 1), the right end of the last burnt bridge being at the point mN. Parameters 
m and A; > assume integer values. 

The dynamics of the system is described by a set of Master equations: 

^^""^^^^^ = wPmN+l,m{t) - {U + P2W) PmN ,m{t) + PlU ^ PmN-l,ni'{t), (29) 

m'=—oo 



for k = i = 0, 

dPmN+kN+{N-l),m{t) 



wPmN+{k+l)N,m{i) + uPmN+kN+(N~2),m(t) 



dt 

+P2Wf{k + l)PmN+(k+l)N,m+k+l{t) - {u + w) PmN+kN+{N-l),m{t) , (30) 

for /c > and i = N — 1 [with f{k + 1) the same as in ([2])], 

^^-^ — — = wPmN+kN+l,m{t) + (1 " PljuPmN+kN -l,m{t) — {u + w)PmN+kN,m{t) , (31) 

for k > 1 and i = 0, and 

■ = wPmN +kN+i+l,m 

[t) + uPmN+kN+i^l,m{t) - {u + w)PmN +kN +i,va 

it) (32) 

for A: > and i = I,- - ■ , - 2. 
We observe that 

RkN+i{t) = ^ P{-m+k)N+i,m{t), {k > 0] 2 = 0, 1 , ■ ■ ■ , - 1) , (33) 

m=— oo 

where RkN+i{t) is the probability for the random walker to be found kN + i sites apart from 
the last burnt bridge at time t, which was used above to find the walker's velocity and is 
given by Eq. ([25]). Plugging ([33D into Eqs. ([29]) - ([32]) results in the equations ([I]) - ([!]) for 
RkN+i, thus obtained by a different method. 

In accordance with p, [lO|, we introduce auxiliary functions SkN+i{t), 

^fc7v+i(t) = 5Z (mAr + /ciV + ^)P^jv+fc7v+i,m(t), {k>Oi = 0,l,--- ,N-1). (34) 

m=—oo 
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The system of equations describing the time evolution of the functions SkN+iit) results from 
Eqs. dMD and Eqs. ([29]) - ([32]). It was obtained that 

dSo{t) 



.^Siit) - {u + p2w)Soit) +piu'^SaN-i{t) +Piu'^RaN-iit) -wRi{t), (35) 

a=l a=l 



^^''^^dt ^^^'^^ ^ wS(^k+l)Nit) + uSkN+{N~2) (t) + P2Wf{k +l)[So- Ro] 

-{U + w)SkN+iN-l){t) - wR(^k+l)N{t) + uRkN+iN-2){t), (36) 



for A; > 0, 

dSkNit) 



wSkN+l(t) + {l-pi)uSkN-l{t) - {u + w)SkN{t) + (1 - Pi)uRkN~l(t) - wRkN+l(t) 



dt 

(37) 

for k > 1 and 

dSkN+z{t) ^ yjS^,^^i_^_^{t)+uSkN+i-i{t)-{u+w)SkN+i{t)+uRkN+i-i{t)-wRkN+i+i{t) (38) 
dt 

for > and i = 1, ■ ■ ■ , - 2. 

At t — >^ oo the solutions of Eqs. (!35|) - (138!) are sought in the form 

S,{t) = a,t + T„ (39) 



where aj and T,- are time-independent coefficients. Plugging Eq. (l39l) into Eqs. ( 1351) - (1381) 
leads to, 

oo 

Wai - (u +P2W^)ao +PlM^aaAr-l = 0, (40) 

«=1 

Wa(fc+i)Ar + MafcAr+(iv_2) + P2W f [k + l)ao - (m + w)akN+{N-l) = (41) 

for k>0, 

wakN+i + (1 - Pi)MafcAr_i - (m + w)akN = (42) 

for k > 1 and 

wafcAT+i+i + ua^Tv+i-i - {u + w)akN+i = (43) 

for > audi = I,-- - ,A^-2. 

Clearly, Eqs. ( HOi) - ( l43l) are identical to the system of equations dH) - (jll) for the functions 
Rj in the t —>■ oo limit, where dRj/dt = 0. Thus their solutions should coincide up to the 
multiplicative constant, namely, 

o-kN+i = CRkN+i, (44) 
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with RkN+i given by Eq. ( l25l) . The normahzation condition ^ ^ RkN+i = 1 imphes that 

fc=0 i=Q 

oo Af-1 

C = ^ ^ akN+i- To find the exphcit expression for C, we utihze the equations for Tj 

k=Q i=0 

obtained by plugging Eq. ([391) into Eqs. fl35l) - fl38ll. 

oo oo 

ao = wTi -{u + P2w)To + piu ^ T^n-i + Piu i?aiv-i - wRi, (45) 

Q=l a=l 

dkN+iN-l) = wT(^k+l)N + uTkN+(N-2) + P2W f {k + l)[ro - Rq] - (m + w;)TfeAr+(Ar_i) 

+Mi?fcAr+(Ar„2) " wR(^k+l)N (46) 

for > 0, 

OfeAT = wTkN+1 + (1 - Pl)uTkN-l - {U + w)TkN + (1 - Pl)uRkN^l - wRkN+1, (47) 

for > 1 and 

OfcAf+j = wTkN+i+1 + uTkN+i-1 — {u + w)TkN+i + uRkN+i-1 — wRkN+i+l, (48) 

00 

for A; > and i = 1, ■ • ■ , — 2. Summing up Eqs. ( l45l) - ( HHI) and using ^ /(/e + 1) = 1 
yields 

oo Af-1 

C = ^ ^ Ofciv+i = w(l - p2)Ro + {u-w) = V, (49) 

fc=0 i=0 

where the walker's velocity V is given by fl28|) . Hence, in accordance with Eq. (jH]) 

O-kN+i = V RkN+i (50) 

where RkN+i and are given by Eqs. fl25]) and fl2S]) . ^ 

Now we are able to obtain the expression for the random walker's velocity [3]. The mean 
position of the particle is given by 

+00 oo A''—! 



= 5Z ^i^^ + kN + i)PmN+kN+i,m{t) = 

m=—OD k=0 i=0 
+00 1 oo Af-1 

^ ^""^ + kN + i)PmN+kN+i,m{t) f = Z Z ^kN+iit), (51) 



fc=0 j=0 l^m=— oo ) k=0 i=0 

which results in the mean velocity 

, oo N-1 , 

d , , , , \ ^ X ^ a 



v = j^{m = T.T.jt^^^^^(t). (52) 

A:=0 i=0 
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In the t ^ oo limit Eqs. (139!) and (H9|) therefore yield 

oo N-1 

y = ^ ^ ttkN+i = w{l - p2)Ro + (u-w) = V. (53) 

fc=0 i=0 

As expected, Eq. ( 153|) reproduces the expression (l28l) for obtained above with the use of 
the reduced chemical kinetic scheme method. 

For the purposes of computing the diffusion coefficient D, functions T/^N+i need to be 
found from (gS]) - dH]) 0, llO]- In <^ - (|48]), akN+i should be expressed according to fISOD . 
In analogy with finding Run+i from ([I]) - (jl]), we start with solving fHHI) . In the expression 
( !25l) for RkN+i, we denote 

^ = 1-/3^ ' (54) 

so that 

i?fciv+. = i?fc7v[l-e(l-/?0], (55) 
where RkN = Rov'', with Rq, y given by ([26]), ([23]). With the use of ([55]) , ([48]) takes the form 

wTkN+i+i + uTkN+i-i-{u + w)TkN+i + {,u-w-V)RkN{l-0 + i^~u-V)RkN^P' = 0. (56) 

We seek the solution of ([56]) in the form 

TkN+i = Ci{k) + C2{k)(3' + iA{k) + iB{k)(3\ (57) 

In ( I57]) . Ci(A;) +C2(/i;)/9* part is the solution of homogeneous equation [the part of ([56]) which 
involves only T^N+i, TkN+i±i], in analogy with Eq. ([8]). Substituting Eq. ([57]) into Eq. ([56]) 
gives these expressions for A{k) and B{k): 

A{k) = ''~''~^ RkN{l-0, B{k:) = ''~'"^^ RuNi. (58) 
u — w u — w 

Plugging Tfcjv, TfcAT+i instead of T^N+i into ([57]) . one can express Ci{k) and C2(/i;) in terms of 
yl(/i;), B{k) as well as T^tv and T^n+i (first two points of the period). Substituting resulting 
expressions for Ci{k) and C2{k) together with the expressions ( [58]) for A{k), B{k) in ([57]) 
gives after some rearrangement 



TkN+i — TkN + dk-. ^ + RkN (1 — 

1 — p u — w 



u-w + V ^ 

-TttkN ^ 

u — w 



1-/3 
1-/3'' 



(59) 
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where 

dk = TkN+i — TkN- (60) 

Eq. fl59|) solves fl56|) [i.e., the equation fj48l) ] for arbitrary T^^r and 9^. Parameters T^n and 
^fc are to be determined from the Eqs. fH6|) and (1471) . To get a less cumbersome result for 
TkN+i, it is convenient to rewrite Eq. fH6l) by adding and subtracting wT^n+n, wRkN+N 
[we note that TkN+N, RkN+N are given by (l59l) and (l55l) with i = N, i.e., TkN+N 7^ 7'(fc+i)Ar, 

RkN+N 7^ -R(fc+l)Ar]: 

afcAr+(Ar_i) = wT(fc+i)iv - wTkN+N + wT^n+N + uTkN+(N-2) + P2Wf{k + l)[To - i?o] 

-(m + M;)TfcAr+(Ar_i) + uRkN+{N-2) - wRkN+N + w{RkN+N - -R(fc+i)Ar), (61) 

which gives 



w(T(fc+i)iv - TfcAT+Ar) + 



i?(,+l)jv) + P2W/(fc + l)[To - i?o] = 0, 



(62) 



as follows from (l48l) which was formally extended to include i = — 1. In fl62|) . the function 
/(A; + 1) is expressed according to Eq. (I7l). Plugging ( !59l) . ( l55l) (with appropriate k, i) into 
(1621) permits us to express 6'^ in terms of T^at in this way, 

1-/? 



{T(k+1)N 



TkN-FRkN+P2y\l-y)To} 



where 



u — w — V 



u — w 



1-0 



1-/3 



u — w + V 



u — w 



1-/3^' 
1-/3 



(63) 



(64) 



It should be mentioned that substituting ( 159|) and ( 155|) into original Eq. ( H6|) leads to a more 
involved expression connecting 6'^ and TkN which turns out to be identical to Eq. (l63l) . as 
was numerically verified. 

To find TfcAT, we substitute fl59|) and fl55l) (with appropriate k, i) in Eq. fHTl) . with 6^ in 
( !59|) replaced according to (!63l) . After some algebra, this results in 



aTfcAT + hTi^k+i)N + cT(fc+2)Ar + {(i + + fy^}Toy'' + xy'' = 0, 



(65) 



where coefficients a, b, c, d, e, / are given by Eqs. ( 1151) - (1201) and 



■ 1-/3 , . 1-/?^-^ 



1-/3^ 



1-/3 



TV 



+ Ro{il-pi)u<P-wy[l-ai-P)] 

+ (l-pi)«[l-e(l-/3^'^)]-V^2/}. (66) 
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In Eq. f l66p it was found that 
u — w — V , 



u — w 







A^- 1 



1-/5 



N-l 



1-/3 



— w + V ^ 
u — w 



iN-l)P 



N-l 



/3- 



N-l 



1-P 



(67) 



and F is given by ([6] 

Comparison between fl65l) and fll4l) shows that they are identical (with replacement R 
T), with the exception of xv'' term in fl65l) . This implies that 



(68) 



with arbitrary constant B solves homogeneous equation [Eq. fl65|) without xy'^ term], since 
y given by (125]) is constructed for R^n = Roy'' to solve specifically Eq. (HM . 
To account for xv'' term, we look for solution of ( l65l) of the form 



r(2) 

-kN 



where 



-h - - 4ac 
2^^ 



(69) 



(70) 



with a, 6, c given by flTSl) - flTTl) : parameter C is to be determined. Parameter 2; is the 
solution of equation cz^ + 6z + a = such that < z < 1. Plugging (169|) into (165|) gives 

C'z'=(a + 6;z + c^^) + (rf + + /y^)^^'^ + XV^ 



0. 



The first term in (ITTi) vanishes since 2; is given by fITU]) . thus (ITTi) yields 



(71) 



(72) 



/nX ~ _____ 

The T^jy given by (169|) with C specified by ( I72|) therefore solves Eq. (165|) . More generally, 
( !65|) is solved by the sum of contributions ( !68l) and ( !69l) . i. e. 



(73) 



where i? remains undetermined. Plugging (1751) in Eq. (jHSD for gives 



/9 



5/(1/ - 1) + C;.'=(;^ - 1) - Fi?oZ/' + PsJ/^l - y)[5 + . (74) 



Substituting (1731) and fl74l) into (159!) results in the final expression for T^n+i, 



TkN+^ = By" + Cz'' + 

+Roy'' 



u — w — V 

u — w 



By\y - 1) + Cz\z - 1) - FR,y'+p,y\l - y)[B + C] 



I — 



/?^ 



/3 



■u — w + 
+ e 



u — w 



ip' - (3- 



(75) 
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Although Eq. (H5|) was not utihzed to obtain (1751) . it was numerically checked that every 
equation in the system fH5|) - fHHl) holds for T^N+i^ RkN+i given by fl75|) and fl25|) . 

To compute the diffusion coefficient D, we need to consider additional auxiliary functions 

+00 

UkN+i{t) = J2 ^"^^ + ''^ + '^^^"-^+kN+i,m{t), k>0, 1 = 0,1, ■■■ ,N -1. (76) 



7Tt= — 00 



A system of equations which determine the time evolution of Uj{t) is derived with the use 
of Eqs. (EED, (El and (1331) for UkN+i{t), SkN+i{t) and RkN+i{t), as well as Eqs. ([29]) - ([321). 
It follows that 

= wUiit) - {u + p2w)Uo{t) - 2wSi{t)+wRi{t)+p^uJ2UaN~i{t) 

a=l 



+2piU^SaN-lit) + PiU^RaN-l{t), (77) 



a=l a=l 



+P2Wf{k + l)Uo{t) - 2wS^k+l)Nit) + 2uSkNHN~2)it) 
-2p2Wf{k + l)S'o(t) + wR(k+l)Nit) + uRkN+(N-2)it) 

+P2wfik + l)Roit) (78) 



for k>0, 

dUkNjt) _ 
dt 

for k > 1, and 

dUkN+i(t) 



wUkN+iit) + (1 - Vi)uUkN-i{t) -{u + w)UkN{t) - 2wSkN+l{t) 

+2(1 - pi)uSkN-i{t) + wRkN+i{t) + (1 - pi)uRkN~i{t) (79) 



wUkN+i+l{t) + uUkN+i-i{t) - (m + w)UkN+i{t) - 2wSkN+i+l{t) 



dt 

+2uSkN+i-i{t) + wRkN+i+i{t) + uRkN+i-i{t) (80) 

for > audi = I,-- - ,A^-2. 

The diffusion constant is to be found from 

D = U\mj^[{x{tf)-{x{t)f]. (81) 
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Using (176!) and summing up Eqs. (177|) - (IHOi) results in 

, , +00 00 N—l 



mN+kN+i,m(t) 
00 TV-l 



m= — 00 A;=0 i=0 

fc=0 i=0 fc=0 i=0 



dUkN+i(t) 



00 7V-1 



+ «;) - w{l - p2)Ro{t) + 2«;(1 - P2)So{t) + 2{u-w)J2Yl 



(82) 



A:=0 i=0 



thus we do not need to find UkN+iif) from the system (1771) - (IHOl) to obtain the diffusion 

00 Af-l 

coefficient. To derive (152]) we used the normahzation conditions, ^ ^ RkN+i{t) = 1 and 

fc=0 i=0 

00 

f{k + 1) = 1. In the t ^ cx) hmit, SkN+i{t) akN+it + TkN+i, Ro{t) Ro, thus ([82]) 

fc=0 

becomes 



00 N-l 



j^{x{tf) = (M+t/7)-w(l-p2)^o+2u»(l-P2)[aot+To]+2(M-M;)EE[afe;v+it+Tfcjv+i]- (83) 



fc=0 i=0 



00 Af-1 

Given that in the stationary-state hmit ^{x{t)) = ^ ^ cikN+i = w{l—p2)Ro + {u — w) = V 

k=0 i=0 

[Eq. (153]) ] and utihzing {x(t)) given by ( 15T]) . it foUows that 



, , 00 Af-l 

-(x(t))^ = 2(x(t))-(x(t)) = 2y(x(t)) = 2V^5^ 5^ W(t - 00) 

k=0 i=0 



00 Af-1 



(84) 



D 



2[w{l - P2)Rq + (m - w)] E E {akN+it + TfcAT+i} . 

k=0 i=0 

Plugging (183]) and (15^ into (IHT]) gives the diffusion constant, 

00 N~l 

u + w) - w{l-p2)Ro + 2w{l -p2)[aot + To] - 2w{l - P2)Ro'^'^{akN+it + TkN+i} 

k=0 i=0 

(85) 

The time-dependent part of (185]) is 

00 N-l 

2w{l - p2)aot - 2w(\ - p2)Ro 

^kN+it I 



1 



fc=0 i=0 



(86) 



00 N-l 

Using ^ CLkN+i = V [Eq. ( 153]) ] and ao = VRo [Eq. (150]) ]. it follows that according to (186]) 
D(t) = 0. Therefore, as anticipated, the diffusion constant does not contain time- dependent 



terms, and it is given by 



2 



00 N-l 



{u + w) - W{1 - P2)Ro + 2w(l - P2)To - 2w{l - P2)Rq E E TkN+i 



k=0 i=0 



(87) 
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In order to get the final expression for diffusion coefficient, it is necessary to calculate 

oo N-1 

E E TkN+i in dHZD. Using Eq. ([75]) for T^N+i yields, 

fc=0 i=0 



oo N-1 



BN CN 



fe=0 1=0 



where 



l-y I- z 



+ 



N 



I -(3^ 1-/3 



-B-C 



RqF 



P2iB + C) 



+A(88) 



A 



Rq \ u — w — V 



\ — y y u — w 



(1-0 



N{N-r 



u — w + V 

+ e 

u — w 



1-/3 



N 



N 



1-/5 
1-/3^ 



(1-/5)2 1-/5 v^' 1-/5 yj i ■ ^^^^ 

In deriving (IHHl) . we used the fact that 0<?/<l,0<2;<l. Next, we consider the last two 
terms in (Ej 



oo N-1 



2w{l - p2)To - 2w{l - p2)Ro J2 Yl = 2^^^ -P2)\B + C-Ro\- 

k=0 i=0 
1 



BN 



y 



CN 



+ 



N 



1-z [l-p^ _ 



-B-C 



RnF 



l-y 



+ P2iB + C) 



+ A 



(90) 



(91) 



where we used (IHHll and Tq = B + C [Eq. (1731) ]. The contribution from terms oc B in (l90ll 
can be shown to be 

' N 1 ' 

l-f3^ ~ l-f3_ 

Utilizing Eq. (l26ll for _Ro, it follows from (19T!) that 5-contribution in (pOl) equals 0, thus 
undetermined constant i? cancels out in flSTl) and it has no effect on the diffusion coefficient. 
Without 5-terms Eq. becomes 



2ti;(l - p2)5 |l - ^ (^iV - (1 - p2)(l - y) 



2^(1-^2^ 



oo N-1 

Tq — Rq TkN+i 
k=0 i=0 

N 1 



1-/3^ 1-/3 



2wil-p2)'!^C-Ro\^- 
RoF 



CN 



l-y 



+ {1-P2)C 



+ A 



(92) 



Plugging fl92|) into fl87|) gives 
1 



D = ^ <! (m + w) - w(l - p2)i?o + 2«;(1 - ^2)6* - 2^(1 - p2)Ro - 



CN 
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In Eq. (p3|l we have P = u/w, and parameters Rq, C, y, z, F, X are given by Eqs. (l26l) . 
( !72|) . ( !23l) . ( iTOl) . ( !64l) . ( !89l) . correspondingly. Other useful parameters [on which D in ( !93|l 
imphcitly depends] are a - /, x, ^, 1^ given by Eqs. (HSl) - ([20]), ([66]), ([M]), and ([281), 
correspondingly. With the help of these parameters Eq. ( [93i) gives the exact expression for 
D as a function of transition rates u, w, probabilities pi, p2, and concentration of bridges 
c = 1/N. 

It was verified numerically for various pi and c values that in the limit of u ^ 1, w ^ 1, 
P2 Eq. (1931) reproduces the diffusion constant obtained in 7] for BBM with u = w = 1 
and p2 = 0. In the limiting case of = 1 we obtained D{u,w,p2, N) in ^ using Derrida 
method [10] and it also agrees with our general result ([93i) in the pi ^ 1 limit, as was 
numerically checked. 



DISCUSSIONS 



To illustrate our findings, we plot the dynamic properties of the molecular motor using 
Eqs. ([28]) and ([93]) . We first consider the case of the unbiased molecular motor with 
transition rates u = w = 1 (Figs. 3-5). Eqs. ([28]) and ([93]) cannot be used directly when 
(3 = u/w = 1. Whereas it is possible to find the u 1, w 1 limit of Eq. ( [28]) (although 
it leads to a cumbersome expression) and thus to plot the velocity, it is problematic in the 
case of Eq. ([93]) . Thus we plotted the diffusion constant for u = 0.999 and w = 1 (see 
Figs. 3(b), 4(b), 5(b)). Using u values closer to 1 generates numerical instability. This is 
a good approximation of the u = w = 1 case, as was judged from comparison with known 
limiting cases. Namely, comparing p2 = case [Figs. 3(b), 5(b)] with the result from |3] 
for u = w = 1 showed the discrepancy in D values of the order of ~ 0.001 for almost 
entire range of parameters c and pi, with the exception of small c ^ 0.02, and pi ^ 0.01, 
where discrepancy exceeded 0.008 and 0.006 correspondingly. For Fig. 4(b), we compared 
Pi = 1 case (not shown) with the corresponding case for u = w = 1 obtained in [9^: typical 
discrepancy in D values between u = 0.999 and u = 1 cases was ~ 0.0005 for all c values 
except c ^ 0.001, where discrepancy exceeded 0.007. 

As anticipated, when the recovery probability p2 ^ I and the presence of bridges has no 
effect, V ^ u — w = and D |(m + w) = 1 for all c (and pi) values (Figs. 3 and 5); the 
same happens in the limit of the burning probability pi ^ (Fig. 4). Increasing pi and 
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Figure 3: Dynamic properties of the unbiased molecular motor with u = w = \ and with the 
probability of burning pi = 0.1: (a) The mean velocity as a function of the concentration of 
bridges for different recovery probabilities; (b) the dispersion as a function of the concentration of 
bridges for different recovery probabilities. 




Figure 4: Dynamic properties of the unbiased molecular motor with u = w = \ and with the 
recovery probability p2 = 0.4: (a) The mean velocity as a function of the concentration of bridges 
for different burning probabilities; (b) the dispersion as a function of the concentration of bridges 
for different burning probabilities. 

the concentration of weak links c leads to increasing velocity [as shown in Figs. 3(a), 4(a), 
5(a)], whereas increasing ■p^ reduces the velocity [Figs. 3(a), 5(a)]. 

The diffusion constant plotted in Figs. 3(b) and 4(b) is a decreasing function of the bridge 
concentration c because the presence of bridges lowers the fluctuations of the motor protein 
(the motor protein cannot cross back the already burned bond). We observed that in the 
limit of low c there is a gap in the dispersion [Figs. 3(b) and 4(b)]: D{c 0) differs from 
the expected D{c — 0) value of \{u + w) = 1. This phenomenon corresponds to a dynamic 
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Figure 5: Dynamic properties of the unbiased molecular motor with u = w = 1 and with the 
concentration of bridges c = 0.2: (a) The mean velocity as a function of the burning probability for 
different recovery probabilities; (b) dispersion as a function of the burning probability for different 
recovery probabilities. 
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Figure 6: Dynamic properties of the backward biased molecular motor with u = 0.3 and w = 0.7, 
and with the burning probability pi = 0.3: (a) The mean velocity as a function of the concentration 
of bridges for different recovery probabilities; (b) the dispersion as a function of the concentration 
of bridges for different recovery probabilities. 

transition between unbiased and biased diffusion regimes as was areued earlier in Ref. 19|. 

n 

Fig. 3(b) is similar to the corresponding plot in Ref. \9] for pi = 1 case, but for pi = 0.1 
the gap is prominent only for small p2 values, and for p2 > 0.1 it practically disappears. We 
note that for P2 = case in Fig. 3(b) the correct c ^ limit must be D(c — 0) = 2/3 
However, it did not reach this value because of the emerging numerical instability for very 
low c ^ 0.001 (see also our discussion above). 

Analysis of Fig. 5(b) shows that as p2 increases, the behavior of diffusion constant 
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Figure 7: Dynamic properties of the backward biased molecular motor with u = 0.3 and w = 0.7, 
and with the recovery probability p2 = 0.1: (a) The mean velocity as a function of the concentration 
of bridges for different burning probabilities; (b) the dispersion as a function of the concentration 
of bridges for different burning probabilities. 




Figure 8: Dynamic properties of the backward biased molecular motor with u = 0.3 and w = 0.7, 
and with the recovery probability p2 = 0.6: (a) The mean velocity as a function of the concentration 
of bridges for different burning probabilities; (b) the dispersion as a function of the concentration 
of bridges for different burning probabilities. 

changes from increasing to decreasing function of pi, with D{pi) developing a minimum for 

P2 ^ 0.1. It should be noted that in 02 = case, D should approach the value of 1/2 for 

n 

Pi according to [7]. We see some discrepancy there which is also due to the numerical 
instability for small pi values. In addition, we observed a gap between D{pi — > 0) = 1/2 for 
P2 = and D{pi 0) = 1 for all nonzero p2 values. 

As another example, we investigated the case of the backward biased motor protein (with 
specific transition rates u = 0.3, w = 0.7), where bridges are inducing the molecular motor 
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Figure 9: Dynamic properties of the backward biased molecular motor with u = 0.3 and w = 0.7, 
and with the concentration of bridges c = 0.2: (a) The mean velocity as a function of the burning 
probability for different recovery probabilities; (b) the dispersion as a function of the burning 
probability for different recovery probabilities. 

to move in the opposite direction (Figs. 6-9). The analysis of dynamic properties shows 
that in the p2 ^ 1 limit (the deterministic bridge recovery) V ^ u ~ w = —0.4 and D — *■ 
!(« + w) = 1/2 for all c and pi values (see Figs. 6 and 9), as it should be. The pi — > limit 
in the u < w case is more complex than in the u > w case (when V ^ u — w, D ^{u + w) 
for all P2 7^ and c values as pi — 0). Namely, ioi u < w there exists p2{u, w) such that for 
non-zero p2 < P2 the dynamic properties V{c) and D{c) exhibit strong c-dependence in the 
Pi — s> limit and they are different from the u — w and ^{u + w) values for non-zero c [Fig. 
7(a), (b)]. Exactly at = (with p2 < P2), V{pi = 0, c) = u — w and D{pi = 0, c) = Mu+w) 



as obtained using the y = 1 root of Eq. fl22|) . given by Eq. fl23l) with "+" rather than "— " 
sign before the square root [the second root of Eq. fl22|) ]. Thus for p2 < P2 there is a 
dynamic transition separating pi = and pi ^ regimes. For u = 0.3, w = 0.7, it was 
found that p2 ~ 0.57. For p2 > P2, V{pi — > 0, c) = V{pi = 0,c) = u — w = —0.4 and 
D{pi — > 0,c) = D{pi = 0,c) = ^{u + w) = 0.5 for all c values [Fig. 8(a), (b)]. Physically this 
implies that bridges with arbitrarily low burning probability strongly affect the dynamics 
of the particle which tends to move in the backward direction provided that the recovery 
probability p2 is less than some critical value; otherwise (in p2 > P2 case) weak links have 
no effect on the particle dynamics as pi — > 0. Thus in the pi limit there is a dynamic 
transition at p2 = P2 separating the regime with p2 < P2 when weak links play a role in the 
motor protein dynamics and the regime where they are irrelevant {p2 > ^2)- It should be 
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noted that the pi ^ hmit in Fig. 9 is in agreement with that in Figs. 7 and 8: p2 ~ 0.57 
plays the same role. The dynamic transition in the pi —>■ limit also takes place for P2 = 
(irreversible burning of weak connections). It separates backward biased pi = regime 
from forward biased regime with small finite pi, when velocity V{c) is positive for all c > 0, 
although V{c) ^ 0, D{c) for all c values as pi — > [see Fig. 9(a),(b) with the specific c 
value]. There are therefore jumps in V{c) and -D(c) at = for all p2 < P2 and all non-zero 
c. 

Increasing the burning probability pi and concentration of weak bonds c reduces the 
magnitude of particle's velocity in the backward direction; the same effect is observed when 
the recovery probability p2 is reduced, as expected [see Figs. 6(a), 7(a), 8(a), 9(a)]. For 
sufficiently large c, pi (and small P2) the velocity V even becomes positive [Figs. 6(a), 7(a)]. 
We observed that for p2 = 0,V is always positive as the burning of weak bonds is irreversible 
in this case [Figs. 6(a), 9(a)]. In that case, V{c = 0) = u — w is different from V{c 0) = 
[Fig. 6(a)]. This effect was also observed in Ref. [9] in the pi = 1 case. 

Diffusion constant demonstrates a more complex behavior, with large fluctuations at 
small c and small p2 which increase with increasing pi [Figs. 6(b) and 7(b)]. Fig. 6(b) with 
Pi = 0.3 is qualitatively similar to the corresponding figure in [9|] with pi = 1, although the 
maxima in D curves are less pronounced than in the Ref. 91]. The shape of the D curve 
differs significantly between Figs. 7(b) and 8(b) when the threshold p2 is crossed. In case 
of P2 = (irreversible bridge burning), fluctuations are reduced (especially at low c) and 
D curve differs substantially from non-zero p2 case [Figs. 6(b) and 9(b)]. As was the case 
with the velocity, for p2 = there is a gap in D separating D{c = 0) = ^{u + w) from 
D{c ^ 0) = [Fig. 6(b)], again illustrating a dynamic transition. For non-zero p2 the 
diffusion constant is D{c ^ 0) = |(m + w), and there is no gap [see Figs. 6(b), 7(b), 8(b)]. 



CONCLUSIONS 



We have presented a comprehensive theoretical method of calculating dynamic properties 
of molecular motors in reversible burnt-bridge models for periodic bridge distribution. It is 
a generalization of the approach developed by us in Ref. 0, Q] for the unbiased molecular 
motors and irreversible burning of bridges. Exact and explicit expressions for mean velocity 
and dispersion have been derived for arbitrary values of parameters u, w, pi, p2 and c. In 
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the known limiting cases ofu = w = l,p2 = and of pi = 1, we have reproduced our earher 
findings g], 3, Q] , thereby confirming the vahdity of our theoretical analysis. Some interesting 
phenomena have been observed as a result of the investigation of dynamic properties of the 
molecular motor in BBM with bridge recovery. It includes dynamic phase transitions and 
reversal of the direction of the motion. In case of the unbiased molecular motor, increasing 
the concentration of bridges c (or lowering the recovery probability P2) with other parameters 
kept fixed results in increasing velocity and decreasing dispersion. However, dependence of 
the dispersion on burning probability pi is more complex; it is determined by the p2 value. 
In the limit of low c, gaps in dispersion plots have been observed for various pi and p2 values, 
indicating the dynamic transition between biased and unbiased regimes. Also, a gap was 
found in the limit of small pi between p2 = and non-zero p2 regimes. Thus our results 
obtained in for pi = 1 with u = w were generalized to cover the full range of pi values. 

For the backward biased molecular motor, increasing c has resulted in slowing down the 
backward movement of the particle, and for sufficiently small p2 (large pi) the direction of 
motion has been even reversed and the velocity became positive. In the limit of small pi, 
a dynamic phase transition separating pi = and — > regimes has been found provided 
that P2 is less than some critical value. For sufficiently small p2, broken bridges infiuence 
the particle's dynamics even if the burning probability pi is infinitesimal. The behavior of 
dispersion as a function of c was non-monotonic for some range of parameters pi and p2, 
with large fiuctuations at small c and small p2- In the case of irreversible bridge burning 
{p2 = 0), we have observed gaps in velocity and dispersion in c — limit (for pi = 0.3), 
with the velocity being positive for all non-zero c values. It suggests that there is a dynamic 
transition at c = separating backward biased and forward biased diffusion. The velocity 
and fluctuations are suppressed for sufficiently small c. Hence our findings in for u < w 
case with pi = 1 have been extended to describe the general case of < < 1. 

The method presented above applies to the case of periodic distribution of weak bonds, 
which is probably realistic for coUagenases bl. As a problem to be addressed in the future 
studies, one can consider BBM with random distribution of bridges [J, 15[ where a different 
theoretical approach must be applied. 
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